Holographic Confinement /Deconfinement Transitions 
in Asymptotically Lifshitz Spacetimes 

Benson Way 



Department of Physics, UCSB, Santa Barbara, CA 93106 
benson@physics.ucsb.edu 



Abstract 

Under a Schcrk-Schwarz compactication, the AdS soliton and blaek branc provide 
a dual description to a confinement /deconfinement phase transition. We extend this 
construction to asymptotically Lifshitz spacetimes. In particular, we show that there 
must be a phase transition between the soliton and black hole if the solutions exist. 
We also construct these solutions numerically and compute the phase diagram. 



1 Introduction 



Since the emergence of the AdS/CFT correspondence [H [21 E], holography has become a 
powerful tool for understanding certain strongly coupled field theories. In particular, there 
are many promising attempts to model interesting condensed matter systems with a gravity 

dual HEIE]. 

Most of these holographic models use a gravity dual which is asymptotically Anti-de Sitter 
(AdS). The field theories described by these models are relativistic and have a conformal 
symmetry in the ultraviolet. However, many condensed matter systems are non-relativistic. 
For instance, quantum critical systems often exhibit a Lifshitz scaling symmetry: 

t \H, X ^ Ax , (1.1) 

where time and space scale anisotropically. Such a scaling symmetry is modeled holograph- 
ically by a Lifshitz metric [7]: 

ds^ = e (-r^'df + ^ + r^dxidA . (1.2) 



This metric exhibits the scaling (1.1) if one also scales r — )■ A~^r. li z = 1, this spacetime is 
AdS in Poincare coordinates with AdS length scale £. 

Attempts were made to understand the holographic dictionary for asymptotically Lifshitz 
spacetimes [Hlin]. To describe these systems at finite temperature, one must construct a black 



hole that asymptotically approaches ( 1.2 ). Indeed, these types of solutions have been studied 

iniioiiiiiiiaiiaiiiiiisiiisiiiz]. 

In this paper, we aim to study a new class of solutions that are asymptotically Lifshitz. 
However, we require a new boundary condition. We simply impose that one of the spatial 
directions be compactified to a circle (with fermions being antiperiodic around this circle, 
which manifestly breaks supersymmetry [IB]). For z = 1, there is a bulk solution with 
this boundary condition given by the AdS soliton 120]. In the case of asymptotically 
AdS^ X solutions which are dual to A/" = 4 super Yang-Mills, this boundary condition 
introduces masses for the fermions and scalars. The AdS soliton is the gravitational dual to 
this state and represents the confining vacuum with a mass gap [T^l [201 121] • This system 
also exhibits a phase transition at higher temperatures into a deconfining phase described by 
a black hole. In the context of holographic superconductors [221 1231 [211 [25] , the AdS soliton 
was used to model an insulator /superconductor transition 
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We would like to generalize the AdS soliton solution to z > 1. We opt to use the massive 
vector model considered in [HI |9] where the holographic dictionary has been developed. 
A similar solution for z = 3 was found in [TB] using a massive gravity theory in three 
bulk dimensions. These solutions will also have the additional advantage that their zero 
temperature limit is well-defined. If z 7^ 1, the Lifshitz spacetime suffers from a curvature 
singularity in the deep infrared in the form of diverging tidal forces [H [71 EH |29] . The zero 
temperature limit of black hole solutions will also suffer from this singularity. While the true 
nature of this singularity has yet to be understood, it may pose a problem for understanding 
the deep IR physics of the theory. In the soliton solutions, the spacetime caps off and there 
is no curvature singularity, even at zero temperature. See [30l l3T] for other ways of resolving 
the singularity. 

In the following section, we present the massive vector model containing Lifshitz solu- 
tions and review the AdS soliton and its confinement /deconfinement transition in z = 1. 
Then in section 3, we present an ansatz which we used to study the equations of motion. 
In section 4, we discuss the thermodynamics of the system where we prove that a confine- 
ment /deconfinement transition exists for any z where solutions also exist. We then construct 
these solutions numerically and compute a phase diagram in section 5. We finish with a few 
closing remarks. 



2 Action and AdS Solutions 

For condensed matter applications, we would like to model a field theory in (2+1) dimensions. 
Because of the extra holographic direction and the compactified circle, we choose to work 
in five bulk dimensions. It is straightforward to generalize our analysis to any number of 
dimensions. Therefore, consider the following action in five dimensions, 

S = Jd'xV^(^R-2A- \f,,F^'' - ^^M^^) , (2.1) 

where F = dA. The equations of motion from this action are 

R,. = ^9^. + If,,F/ + "^A.A, - ^F,^F^-g,, 
Vf.F^'' = m^A" . (2.2) 
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The action (2.1) is the five-dimensional version of the one considered in [HI E], where much 
of the holographic dictionary has been developed. 



One solution to (2.2) is the Lifshitz spacetime 



(2.3) 



This metric is equivalent to the one presented in (1.2) with r = 1/p. The other fields and 
constants are given by 



A 



A 



m 



dt 



\ z 

z"^ + 2z + 9 



2P 



?>z 



(2.4) 



Now let us set 2; = 1 in (2.4). Then the equations of motion reduce to those of pure 



Einstein gravity with a negative cosmological constant. We therefore have the usual asymp- 
totically AdS solutions. In particular, we have the planar black hole, or black brane metric 



ds' 



dn^ 

f{p)dt^ + + drj^ + dx^ + dy' 



fip) = 1 - ^ , (2.5) 
P+ 



\ " ' fip) 

where p+ is the horizon radius. This black hole has temperature and free energy densitjj^ 



7rp+ 



P+ 



(2.6) 



We can also make the coordinate r] in (2.5) periodic with any period we wish. If we then 



perform a double Wick rotation {dt — )■ idr], di] — )■ idt), we obtain the AdS soliton [121 [20] 

P 



ds' 



P 



fip) 



(2.7) 



In order to avoid a conical singularity at p = the coordinate r]'m. (2.7) must have a period 



7 = Trpo 



(2. 



^Here, we are using units where IGttG — 1. 
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The AdS soliton takes the shape of a cigar with a tip at p = Pq. Since the spacetime only 
exists for p < po, the dual field theory is confining and has a mass gap |19]. We can consider 
the AdS soliton at any temperature with a free energy density given by 

Together, the black hole (with the coordinate rj having period 7) and the soliton are 
geometries with the same asymptotics (M^'^ x S*^). At finite temperature, they form two 
competing phases. These solutions have the same Euclidean geometries, differing only in the 
interpretation of which is Euclidean time and which is a spatial direction. To find the 



critical temperature, we compare the free energies (2.6) and (2.9). Since the free energies 
are computed from the Euclidean action, we find, unsurprisingly, a phase transition when 
T = 1/7 with the soliton phase being preferred at lower temperature. This is a first-order 
phase transition, analogous to the Hawking-Page transitiorj^ [32]. In the dual field theory, 
this is a confinement /deconfinement transition. 

We would now like to extend these spacetimes and the corresponding phase transition 
to asymptotically Lifshitz solutions with z > 1. Unlike the z = 1 case (pure AdS), the pure 



Lifshitz metric (2.3) is not invariant under a double Wick rotation. In parciular, the analytic 
continuation of a black hole solution would give a spacetime with different asymptotics. 
Therefore, the black hole solutions and solitons must be constructed independentl}]^ 



3 Ansatz and Boundary Conditions 

Black hole solutions in this model have been constructed before in p ^[TTl[T2] . but we include 
the analysis here for completeness. To find asymptotically Lifshitz black hole and soliton 

^An analogous transition between asymptotically Lifshitz black holes and thermal Lifshitz was considered 
in [H] 

•^A double Wick rotation of a Lifshitz black hole would yield an anisotropic soliton. Similarly, a double 
Wick rotation of a Lifshitz soliton will yield an anisotropic black hole. While these solutions may be of 
interest, such a transformation would make the vector field in our model complex, so we do not discuss them 
here. 
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solutions, consider the following ansatz: 



^^'-'^F{p)dt' + ^ + H{p)dr]^ + dx' + dy^ 



-P 



Rip) 



A 



e 2{z-i) 



(3.1) 



with A and m as given in (2.4). For now, we are only concerned with analytic expressions. 



For numerics, we will need to modify this ansatz, but we defer the discussion of performing 
numerics to section 5. 

The equations of motion can be reduced to 



F" + 



H" + 



+ 



F' R' H' A{z-l) 
p 



F' R' H' 2(z - 1 

1 ^ ^ 1 

F R H p 



F' + 



F' R' H' 2(z + 2) 

1 1 ^ L 

F R H p 



H' + 



z + 2 [ R' H' 2z{z - 1) 
p [li ^ll^ p{z + 2) 

2(z-l)(p0--#)2 
zp"^ 

F' R! 2(^-1)02 2(2-1)" 

y ~ ^ ^ Jrf ~p _ 

'F' H' R' Q r 1 

l^~ll^^^py'~R 
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P 

z 

2p 



= 

H = 




R 



2zH[{z^ + 2z + 9)F + 3{z - 1)0^] 



zp{pH' - QH)F' + 2[- z{z + 2)pFH' + {Qz{z + 1)F + (z - l)(p0' - #)2)/f] 



(3.2) 



Note that the function R can be expressed in terms of the other functions. To preserve 
planar symmetry in the black holes, we set i/ = 1. There are thus two second-order ODEs 
in F and 6 for the black hole and three second-order ODEs in F, H, and 6 for the soliton. 



Our ansatz (3.1 ) gives a one dimensional reduced Lagrangian that generates the equations 



of motion (3.2). As in [121 E]) we can study the symmetries of this Lagrangian to find the 
following conserved quantity: 



d 



and so 



zpFH' + 3 [ - zpF' + 2{z - l)p00' + 2{z - l)z{F - 0^)] h 



2 + 3 



3zp 



zpFH' + 3 [ - zpF' + 2{z - l)p00' + 2{z - l)z{F 



HF 



(3.3) 



3zp 



2+3 



H Ij^ 
HF 



Co, (3.4) 
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for some constant Cq. Of course, one can also verify that this quantity is conserved via the 



equations of motion (3.2). This conserved quantity will be a key component in our derivation 
of a Smarr relation in section 4. 

Let us now discuss the boundary conditions, beginning with the boundary conditions at 
p = p+ or p = corresponding respectively to the horizon and the of the tip of the soliton. 
A black hole solution requires that F (and R) vanish linearly as p — )■ p+. Regularity also 
requires that vanish linearly. A soliton requires that H (and R) vanish linearly as p — )■ po. 
Regularity then imposes conditions on F'{pq) and 0'(po)- 

Now we turn to boundary conditions near p — )■ 0. We demand that the functions ap- 



proach the Lifshitz metric (2.3) as p — t- 0. Then for small p, any deviation from Lifshitz is 



perturbative. Let us now linearize the equations about the Lifshitz solution 

F(p) = l + eFi(p) 
Hip) = l + eH,ip) 

0(p) = l + e0i(p), (3.5) 

where Hi = for the black hole. From the linearized equations of motion, we find that the 
functions have the following behavior near p — t- 0: 

Hi ~ bp'+^ 

01 



where X± = z + 3 ± ^/9z'^ — 262; + 33. See figure [I] for a plot of these powers as a function of 
z. The c coefficients are related to the others via 

_ 3a{2z^ + z + 9) - bz{z + 2){z + 3) _ za±X± 

6(2-3)(^-l) ' ~ ~A{z-3){z-l) ■ ^ ' 

Following the analysis in [S], we also require that the p^^^"^ terms vanish. For z > 3, these 
terms diverg^ as p — )• 0; for 1 < z < 3, these terms can be interpreted as boundary data 
for the vector field [8] . As we will see in the next section, the terms proportional to p^^^ are 
related to the energy density. 

For any given z, there are thus a family of black holes and solitons parametrized by p+ 



*The 2 = 3 case introduces logarithms. 
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Figure 1: A plot of the various powers p appear in (3.6). The dotted red hne is A_/2, the blue 
dashed line is A+/2, and the solid black line is z + 3. 



and po; respectively. These solutions are related to each other by a useful scaling symmetry: 



Given a black hole solution with horizon at p = p+, the above transformation will give 
another solution with horizon at A^^p+. Similarly, a soliton with tip at p = po can be scaled 
to give a soliton with tip at p = A^^po- 

4 Thermodynamics 

In this section, we discuss the thermodynamics of the black holes and solitons. In particular, 
we will derive the thermodynamic Smarr relation found in pEfj and a similar relation for 
the soliton. This will allow us to prove the existence of a phase transition, so long as both 
solutions exist. For simplicity, we will henceforth set IGtcG = i = 1. 

We begin by computing the energy density. The expression for the energy density in 
asymptotically Lifshitz spacetimes has been worked out in four dimensions in [S] from a 
non-relativistic boundary stress tensor complex. Extending their work to five dimensions. 



t^XH 



{p,x,y,r]) -> X{p,x,y,r]) . 
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the energy density is given by the p — )■ hmit of the expression: 



£ = 2s\ - s'At = -2p^'su + p 



2(z-l] 



-St , 



(4.1) 



where 



-h 



2{z-\\ 



-A^A^Y^I'^A, 



(4.2) 



Here, h^J^y is an induced metric at some constant p; is a unit vector normal to the boundary 
and directed outwards; and vr^j, = K^^ — Kh^^^, where i^^j, = V(^n^) is the extrinsic curvature 
of the boundary. 



Given the behavior near the boundary (3.5) and (3.6), the energy density is 



S= ^±l(^_sa + zb) . 



(4.3) 



We would like to relate this quantity to other quantities at the horizon or tip. As in [14|, we 



do this by evaluating the conserved quantity Co given by (3.4). At the boundary, we find 
that Co satisfies 



£ 



-Co. 



(4.4) 



3 + z 

Fo derive the Smarr relation for the black holes, we evaluate Cq at the horizon and then 



relate it to (4.4). At the horizon. 



Co 



2 + 2 



p 



F\p+)R\p+) = Ts , 



(4.5) 



where T is the black hole temperature, and s is the entropy density. Together with (4.4) 
this implies the Smarr relation 



8^ 



bh 



(4.6) 



^ + 3' 

which is the five dimensional version of the same thermodynamic relation derived previously 
in [H]. 
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We can repeat this procedure for the sohton. At the tip, 



^F{p,)R{p,)H'{po) . 



(4.7) 



By analogy with (4.5), we can rewrite (4.7) as Cq = — r/37, where 7 is the period of the 



sohton, and r is the tension. If we were to double Wick rotate the soliton, r would be the 
entropy density of the resulting anisotropic black hole, and I/7 would be its temperature. 



From this and (4.4), we find that 



-2 + 87 

We can compute the free energy density via 

J^ = S-Ts . 



(4. 



(4.9) 



Since the soliton has vanishing entropy, we see immediately from (4.8) that the free energy 



is negative, and so the soliton is always preferred over thermal Lifshitz, as expected. The 
free energy of the black holes are also negative and given by J-" = —zTs/{z + 3). Equating 
the two free energies, we see that there is a phase transition when 



Ts 



T 

^7 



(4.10) 



From the scaling relations (3.8), one can show that decreasing the black hole temperature 



also decreases the entropy density, so there always exists a temperature where (4.10) can be 



satisfied. We therefore conclude that if soliton and black hole solutions exist for a given z, 
there must always be a phase transition between them. 



5 Numerics and Results 

Now we construct the black holes and solitons numerically. To improve numerics, we must 



make modifications to our ansatz (3.1). We must satisfy all of the boundary conditions and 
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keep the functions at least twice differentiable. For the black holes, we set 



Hip) 



1 



0(p) = (l-/^+')(l + P%(p)) 



(5.1) 



The horizon is located at p = p+ = 1. The scaling (3.8) lets us obtain solutions for other 
values of p+. Here, the equations of motion give boundary conditions on f'p{l) and /^(l). 
At the boundary, we require /f(0) = /^(O) = 0. This ensures that the A_ terms in (3.6) 



vanish. The power of 2^; + 2 was chosen so that the other modes in (3.6) are unaffected. 
For the solitons, we set 



F{p) = l-p^gp{p) 
Hip) = 1 - p'^'gnip) 
0(p) = l-p^ 



(5.2) 



If we set 5'_f/(l) = 1, the tip is located at p = 1. The equations of motion then give boundary 
conditions on g'p{l) and g'^{l)- As in the black holes, we require (7^(0) = gni^) = ^^^(O) = 
at the boundary. 



Based on the behavior near p — i- given by (3.6), all of our functions are guaranteed to 
be at least twice differentiable. 

We opt to use a Newton- Raphson relaxation procedure on a finite difference grid. Because 
of the inherent non- analytic behavior of the functions, particularly at p = 0, we use a second- 
order central differencing in the interior of the grid and forward and backward differencing 
on the edges, rather than some higher order method. 

As a test of convergence of our code, we use the quantities 



|1 - iNifp, f<j>)/lN+l{fF-, f4>)\ 



|1 - lN{gF,gH,g4>)/lN+i{gF,gH,gcf>)\ 



(5.3) 



where iNifi) is the average of the integrals of the functions /j, computed using a grid of 
points. As we see in the left plot in figure [2| we find the expected quadratic convergence. 



We can also use the relations (4.6) and (4.8) to test the accuracy of our numerics. We 



compute the energy density from the boundary Sq using (4.3). We then compute the energy 



from quantities at the horizon or tip Sq using (4.6) and (4.8). A measure of our accuracy is 



then given by Max(|l — £q/£o\,\1 — £o/£q\). The result of this error is plotted in figure |2] on 
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Figure 2: Left: test of quadratic convergence for the black hole by computing (5.3) at z = 1.3. The 



soliton is similar with a fit of —0.107137— 1.99508x. Right: the error in the Smarr relation 



(4.6) and the relation (4.8) in red and blue, respectively. Both plots are computed with 



z = 1.3. 



the right. We find that these relations are satisfied to 0.1% accuracy. 

Now we compute the phase diagram using the following procedure. First, we must fix a 



scale for the period of the spatial circle. We use (3.8) to set the period of all of our soliton 
solutions to 7 = TT and compute the free energy of these solutions. For any given z, we then 
compute the free energy of the black hole and scale the solution so that it has the same free 
energy as the corresponding soliton. Finally, we compute the temperature of the resulting 
black hole, giving us the critical temperature for the phase transition. The result is plotted 
in figure [3] We note that as z — )■ 1, the critical temperature approaches the expected value 
of l/vr. 



6 Discussion 

Using a massive vector model in five dimensions, we have shown that a confinement /deconfinement 
transition exists whenever black hole and soliton solutions exist. It would be interesting to 
see if this analysis can be extended to other models that include Lifshitz metrics, especially 
models with embeddings in supergravity or string theory. 

We have also constructed Lifshitz solitons and black brane solutions numerically and 
computed their phase diagram. It seems at first glance that the critical temperature is 
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1.0 1.5 2.0 2.5 3.0 3.5 4.0 



z 

Figure 3: The phase diagram for the Lifshitz black hole/sohton system. The period of the sohton 
has been scaled to 7 = vr, and we have set 16nG = i = 1. 

approaching zero as 2; — )■ 00. It is unclear whether the temperature continues to decrease for 
larger z, especially since numerics become more difficult (as can be seen from the right plot 
of figure [2]). Perhaps something can be said from the 2; — )■ 00 limit of Lifshitz, AdS2 x M^^^. 
However, it is not clear how such a limit should translate to our black holes or solitons. 

One can consider holographic superconducting models of the Lifshitz solitons. This may 
lead to a superconductor/insulator transition with z > 1 Lifshitz scaling. In the z = 1 
case, the AdS soliton represents an insulator, and can be made into a superconductor by 
varying a chemical potential Together with the confinement /deconfinement transition, 
the system has a rich phase structure [27]. One might suspect that these properties can be 
extended to 2; > 1. 
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